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Abstract 


Those terms in the thermal conduction equation which arise 
due to infrared emissions and absorptions by atomic oxygen and 
carbon monoxide are investigated. The purpose of this investi- 
gation is to develop general expressions for the net emission 
by 0 and CO which take account of absorption from the planetary 
black-body spectrum and from other regions of the thermosphere 
and are thus valid at all optical depths in the thermosphere. 

The radiative terms are developed under the assumptions 
of local thermodynamic equilibrium and an Eddington approxima- 
tion for the radiation field. The results are compared to the 
Bates' approximation to the radiative loss terms for a Martian 
model atmosphere which is optically thick in both 0 and CO. It 
is found that the net emission differs from the Bates' expres- 
sions by more than an order of magnitude near unit optical 
depth and that absorption causes the emission to be somewhat 
less than the Bates' approximation in optically thin regions 
of the high thermosphere. 

The expressions developed for the radiative terms permit 
the evaluation of net heating as well as cooling in the O and 
CO emission lines. For the model thermosphere used to compare 
the net emission with the Bates' cooling terms it was found 


that near the mesopause there was a heat source of about 
4/3 

10 ev/cm sec due to heating in the 0 62u line and a lesser 
source of about 10-100 ev/cm sec due to heating in the CO 
rotational lines. 
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The recent advent of direct observation of planetary 
atmospheres by fly-by space probes has given great impetus to 
theoretical studies of the upper atmospheres of the terrestrial 
planets. The occulta tion experiment carried out by the Mariner 
IV spacecraft on the Martian atmosphere has generated a great 
deal of discussion concerning the physical conditions which 
exist in that atmosphere. Since the data obtained fall far 
short of strictly delimiting these conditions, there has been 
wide latitude for theoretical speculations. The observed ioni- 
zation profile has variously been interpreted as an E, F^, or 
F 2 layer, and calculated exospheric temperatures have been 
published ranging from 130° to 400°K (Fjeldbo et al. , 1966; 
McElroy, 1967; Chamberlain and McElroy, 1966) . 

An important aspect of the theoretical computation of the 
thermospheric temperature profile is the effect of radiative 
cooling in the infrared by constituents present in the thermo- 
sphere. The most important radiators in the atmospheres of the 
terrestrial planets are C0 2 , CO, and 0. Radiative losses from 
CO will generally occur near the bottom-of the thermosphere 
below the altitude of solar ultraviolet heating,, whereas CO 
and O may radiate effectively at higher altitudes depending, 
of course, on the amount of these constituents present. 


The usual practice in studying the effect of these last 
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two constituents in cooling a thermosphere has been to utilize 
an approximation introduced by Bates (1951) . This approxima- 
tion takes unit optical depth in a constituent as a dividing 
line above which every emitted photon escapes the thermosphere 
and is thus effective in cooling and below which every emitted 
photon is reabsorbed locally and is thus ineffective in cooling. 
Such an approximation distorts the functional form of the loss 
terms in the thermal conduction equation with consequent error 
in the computed temperature profile. The magnitude of this 
error would, however, be difficult to estimate qualitatively. 

In addition, the Bates' approximation neglects the possibility 
of heating in the radiatively active constituents. This heating 
can occur if the effective temperature of the planet is roughly 
twice the local temperature at the altitude where 0 and CO will 
be formed via photodissociation of CO^ and O^. 

Published models of the Martian atmosphere indicate tem- 
peratures from about 130°K (McElroy, 1967) to as low as 60°K 
(Pieldbo et al. , 1966) in the altitude region where photodis- 
sociation would occur. Since the planetary effective tempera- 
ture is 230°K, the possible effect of heating in the radiative 
constituents must be considered. 

The object of this paper is to develop expressions for the 
radiative losses by 0 and CO which take account of absorption 
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in the emission lines. The use of such expressions in thermo- 
spheric energy balance studies should alleviate doubts which 
may arise concerning the accuracy of such calculations when 
the Bates' approximation is employed in optically thick thermo- 
spheres or in models having a run of temperatures in the lower 
thermosphere which are less than the planetary effective 
temperature. 

The Thermal Conduction Equation 

The most important processes determining the thermospheric 
temperature are the deposition of solar ultraviolet energy and 
the loss of this energy via thermal conduction and infrared 
radiation. These processes are described by the thermal con- 
duction equation which we will write in the following one- 
dimensional, time-independent form: 

i(?«T) ji) + X h, it F.) = - pfe) u> 

i 

h(T) is the thermal conductivity, and Q(B).the ultraviolet 
* 

heat source. The net emission of radiation in a specified 
transition is given by the divergence of the photon flux 
multiplied by the energy of the transition hv. The summation 
in (1) is taken over the possible radiative transitions. The 
factor tt in the emission term arises from the definition of 
the flux as (Chandrasekhar, 1960) 




ip J^L d-fL 


where I is the photon intensity, w the cosine of the polar 
angle, and dn an element of solid angle. The integral on the 
right hand side of (2) gives the number of photons/cm 2 sec 
crossing an area in the radiation beam, and it is the diver- 
gence of this quantity which is of interest. 


Derivation of the Radiative Terms 

The divergence of the radiative flux may be obtained in a 

straightforward way from radiative transfer considerations 

once the appropriate transfer equation has been found. To 

derive this equation we consider a system consisting of two 

energy levels and assume that the spectral line arising from 

transitions between these levels is Doppler broadened. Such 

a system is appropriate to both the ground state • triplet of 

atomic oxygen and to the rotational states of carbon monoxide. 

There are actually two possible transitions within the oxygen 

3 3 

triplet, but the 0( P^) -♦ 0( P^) will be much more important 
3 3 

than 0 ( Pg) -* 0( P^) , since its Einstein coefficient for 

spontaneous emission is about 5 times greater, the statistical 

3 3 / 

weights of and P^ are 1 and 3 respectively, and the 

3 

population of the higher-energy P^ state will be smaller than 
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3 

that of P^. There are a number of excited rotational levels 
for the CO molecule, but only the dipole transitions between 
adjacent states will be of importance; hence the summation 
indicated in (1) may be taken over the rotational states once . 
the flux divergence arising from transitions between adjacent 
states has been found. The radiative emission under discussion 
occurs in planetary thermospheres under conditions of very 
low pressure and the spectral lines are Doppler broadened. In 
what follows we denote the upper level of the two- level system 
by a subscript 2 and the lower level by a subscript 1. 

The number of molecules which undergo transitions 2 -* 1 
and emit a photon in the frequency interval dv will be denoted 
by ^^ 2 , where N 2 is the upper state population. The net 
decrease in the upper state population in a time interval dt 
is -d(6^N 2 ). Equating this to the corresponding number of 
photons added to the radiation field we have 

. - -d(Uk)' 

where dl^ is the increase in photon intensity in the frequency 
interval dv, da is an elementary area across which the radiation 

flows, and the integration is over solid angle. If we multiply 

J 

and divide by an element of path length in the medium, ds, on 
the right hand side of (3) and integrate over frequency we 


Jsl (XL, Jcr dt Jvj 
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obtain 

~ ~lt jjjl j<r j L <lv 

0 

which relates the rate of change of number density in the upper 
state to the change in intensity of the radiation field. For 
a two-level system we may also write 


dt 


- fix Ali + ^ P - B a ^ 


(5) 


where p^ is the density of the radiation field traversing the 
medium, and are the Einstein coefficients for spon- 
taneous and induced emission, is the Einstein absorption 

coefficient, and n^ and n are the number densities in the 
lower and upper energy levels. Since we have assumed Doppler 
broadening, the intensity may be written 


1=1 e 




(6) 


where I is the intensity at the center of the line (at 

v 0 

frequency v ), and a_ = (v /c) (2kT/m) 2 is the Doppler width of 
O D O 

% 

the line. 

We now introduce the assumption of local thermodynamic 
equilibrium (L.T.E.) . • If the radiation field and the gas are 
nearly in equilibrium, the population of the various energy 
states are related through the Boltzmann distribution. In 
particular we have 


J 
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n, _ % 


n, u C -■ ' 't " £fr) 

where g, and g are the statistical weights of the corresponging 

X a 

states, e_. is the energy difference, e - e., between states 
2 x Z x 

2 and 1, e 0 is the energy (above the ground state) of state 2, 
n is the total number density of the gas, and Z(T) is the parti- 
tion function for the system under consideration. The relations 
between the Einstein coefficients are 

a 


- Vi i 


tl. = 


V „ 


e 


n 


(7) 


J*_ _ _c_(c_ 


Jk. = Ja. 

; e» % 


( 8 ) 


Finally, the radiation density and intensity are related by 




(9) 


where I is the mean intensity of the radiation field. We now 
write the element of path length in the medium as ds = uda, 
where da defines the vertical direction. Equating the right 
hand side of (4) and (5) , performing the integration over 

ft 

frequencies using (6) , and making use of (2) , (7) , (8) and (9) , 

we obtain, after some algebraic manipulation, 


1 Jk _ 
i Jt 


= - pj n + s 0 n 


(10) 


where 


) 
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i_ JJa (_m We i ^ T t 
D 2 ( 7 ) (^r 

is the Doppler absorption cross section in the line, and 


( 11 ) 


r _ _jL_ Mil f 01 fjL ) | so-^/kT 

D Z(7) \2kl J \ P J 


(12) 


is the emission coefficient. The ratio e D /a D is the Plank 
distribution (for photons), as it must be since L.T.E. has 
been assumed. The frequency appearing in (10) - (12) is the 
frequency at the center of the line. Equation (10) is the 
standard form of the first moment of the radiative transfer 
equation. The solution of the radiative transfer equation 
for the mean intensity 1^ and flux are discussed in standard 
works such as Chandrasekhar (1960) or Kourganoff (1963) and 
need not be re-derived here. If we assume the adequacy of the 
Eddington approximation, i.e., that the radiation field may be 
divided into upward and downward moving streams of radiation 
having different intensities, then the solution for the mean 


intensity is 

I Jr) = i 1 EJz-r) + i 


,r 0 


B(+) E (w) 


where r is the optical depth in the line defined by 
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JT = - tt 07 <Ji 


(14) 


and E n and E are special cases of the integro-exponential 


Junction defined by 


E„M - IV 


The total optical thickness of the atmosphere in the spectral 
line under consideration is denoted by t . I Q represents the 
flux incident from below, i.e., from the planetary spectrum. 


B(t) appearing under the integral in (13) is the source func- 
tion, which is the Planck function under the L.T.E . assumption. 
It is worth noting that in the thermosphere conditions depart 


radically from radiative equilibrium, and we cannot assume 
I = b, which would make the flux divergence zero. 

The energy from the planetary spectrum which interacts 
locally is that contained in a Doppler width of the spectral 


line. The mean intensity I Q is thus 


i. - a rj c-^ rM(T) 


where P (T ) is the Planck distribution characteristic of the 
V E 

planet, T._, being the planetary effective temperature. If we 
take as approximately constant over the Doppler width of 


the line , then 
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L ^ m ^(Z) *o(r) 


where T is the local temperature at the altitude where I v ( t ) 
is being evaluated. The emission represented by the Plank 
function in the integrand of (13) has a form similar to I q 
except that T_, is replaced by T(t), the local temperature at 
the emitting altitude. 

We may now evaluate the flux divergence for the 0 62 y 
line and the CO rotational lines. Substituting the expression 
for I v (t) of (13) into (10), using (17) for I Q and a similar 
expression for B(t), and evaluating cr D of (11) for the ground 
state triplet of atomic oxygen, we find the divergence of the 
energy flux to be 

The coefficient of the term in brackets is 

-tp’W/T 

ftoxU) s" 1 3e' iif/r -^ e 

which is just the radiative loss term given by the Bates' 

1 2 

approximation. The terms A^(T , T) and A^(T) are 




I £1 


_ | 


L ^ - 1 ) 


( 20 ) 
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^v/ne) | 

A 0 * ( T ) = Tj qUv/TU) _ I t,(l t-V(r)l)dt (21) 

Tries e absorption terms will be discussed after an expression 
similar to (18) is obtained for the CO rotational lines. 

For the CO lines we are interested in the divergence of 
the sum of the fluxes arising from all transitions between 
adjacent rotational states J -♦ J-l, where J is the rotational 
angular momentum quantum number of the upper state. The 
Einstein coefficient for spontaneous emission between these 
states is (based on a rigid-rotator model) 


/It t-, = <T TrT 


I 


(22) 


where d is the electric dipole moment for CO, and B the rota- 
tional constant given by B ="h/4ncl, where I is the moment of 


inertia of the molecule and c the velocity of light. The 
partition function for the rotational states is 


Z(T) = L UTr))e~^r 

X 


r(Ti-l) 



( 23 ) 


The approximate equality in (23) is obtained by replacing the 
sum over J states by an integral. 0 is the characteristic 

3T 

\ 

rotational temperature for CO 


The replacement of the summation over J by an integration is 
valid for T » 0 . Using (22) and (23) we find that the 
expression for the Doppler absorption cross section in the 
J - J-l line is 



M? -ill ( m Jh 

3 r, (2V17 ~ 


J 


C~ T T 


0 



(25) 


We may now follow a procedure similar to that used to find 
the flux divergence in the 62 u O line to obtain the flux 
divergence for a J -* J-l rotational transition in CO. The 
expression so obtained is then summed over J , but the summa- 
tion is replaced by an integration as in the evaluation of the 

partition function, and terms in J are neglected relative to 
2 

terms in J . The result we obtain is 


• ^ = -floM 

J-i J 


v/hore 

Uij). qtcffMdl 

1 

s' 

- 2.5s-io- ,s n,A2)r 


(26) 


3 

is the radiative loss (in ergs/cm sec) given in the Bates' 


( 27 ) 
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approximation, and 



(23) 




[l(t)/i (r)j E,(lt-r60!)c/t 


(29) 


The optical depth, t, is taken as the optical depth in the 
strongest line. The J dependence of the optical depth has not 
been included in the J-sums. 


Comparison with the Bates ' Approximation 

Bates (1951) wrote the rate of radiative loss in the 62 u 0 

line as 

ft = eA n ( l - Vp) < 30 > 

where s is the energy emitted in a transition, A the Einstein 
coefficient for spontaneous emission, and n the number density 
of the emitting state, which is readily related to the total 
species number density through the partition function, ct/p 
is the ratio of the energy of the radiation field which lies 
within the Doppler width of the line to the energy which would 
lie within this width if the radiation field and gas were in 
thermodynamic equilibrium. The expressions derived above for 
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the nat emission in the O, (18), and CO, (26), lines are of 
this form. a/? is the sum of the two terms A 1 and A 2 in these 

0'Z>‘ tii. t 

The first term in the brackets of (18) and (26), a\ 

represents absorption in the line or lines from the planetary 

spactrum.. As noted by Bates, near t = t q this term is roughly 

T_,/2T, where T is the local temperature at t . This follows 
r, o 

directly from (28) for the CO lines and holds approximately 

for the 62^ 0 line if the local temperature is greater than 
o 

228 K, as can be seen by expanding the exponential terms in 
(20). Thus, if the planetary effective temperature is more 
than twice the temperature at the altitude where 0 or CO first 
appears, there can be a net heating due to this term alone. 

Tne second term in (18) and (26), A , represents absorp- 
tion at a given level of the radiation emitted at all other 
levels of the atmosphere, the radiation from each emitting 
layer being attenuated along the optical path between these 
levels . 

If the atmosphere is optically thick, the term represent- 
ing absorption from the planetary spectrum rapidly approaches 
zero at short distances above the CO^ dissociation level. There 
may thus be a large peak of heating due to this absorption near 
the dissociation level, but none of the energy from the planetary 


17 


spectrum will penetrate far into the thermosphere. The term 
representing absorption at a given level due to radiation 
from all other levels does not vanish even at high altitudes. 
If we denote by T the temperature high in the thermosphere 

X 

where r 0, and by a moan temperature at the levels from 
which absorption occurs, then this second term is of the order 


TT ih(T„ 


where r is a mean optical depth from which absorption occurs. 

rx 

In practice, for optically thick atmospheres, we may have 

T /T ^ 1, and t is surely < 1. Thus, this absorption term 
K X M 

should attain some value between, roughly, .08 and .50 at high 

altitudes. The first figure follows from setting t = .7, the 

second from t. = 0. The radiative losses given by (18) and (26) 
M 

therefore do not approach the Bates ' approximation exactly at 
high altitudes, but are somewhat less due to the effects of 
reabsorpti on . 

If the atmosphere is optically thin in the radiating 
constituents, the effect of absorption from the planetary 
spectrum will be felt higher in the thermosphere and thus 
constitutes a heat source which should be included in thermal 


:alculations . The second absorption term will be less 
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important at high altitudes in this case since T , will refer 
“ M 

to a temperature nearer the mesopause, and we will have T^T < 1. 

Figure 1 shows a comparison of the Bates' radiative loss 

term with the loss given by (18) and (26) . These are computed 

qa uVia basis of a modal atmosphere for the planet Mars. This 

is an optically thick model in both 0 and CO, unit optical 

depths in these constituents occurring at 121 and 128 km, respec- 

11 -3 

tively. The model assumes densities of 7 x 10 cm of both 
O and CO at 90 km, which results in optical thicknesses of 
about 6 and 49, respectively. The temperature profile above 
this point is taken as 

T(?) = T n *( T.-T„) 

where T , the temperature at 90 km, is taken as 125°K, and an 
m 

exospheric temperature T = 300°I< has been assumed. The con- 
stent C may be evaluated by assuming a T(z) value at any 
altitude above 90 km. For the computations represented in 
Figure 1, we have taken T(120) = 160°K. The integrals which 
appear in (21) and (29) have been evaluated with a two-point 
quadrature formula, the weights and divisions being computed 
for each value of t. The method is discussed by Chandrasekhar 
(1960) . For the model in Figure 1 we see that there is a 
heating in the CO lines which peaks at 110 km and has a 


-c 


-U-tJ/c 


(32) 
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2 

magnitude of about 32 ev/cai . This is due entirely to absorp- 
tion from the surrounding atmosphere, i.e., to the term A 2 in 
( 25 ) . Absorption from the planetary infrared spectrum occurs 
within the first kilometer above 90 km and is indicated, with 
exaggerated width, at the extreme left hand side of the figure. 

The heat source provided by this term, A 1 , is of the same 
2 

magnitude as A . Both of these heating terras are no doubt of 
trivial importance insofar as determination of the temperature 
in the 90-120 km region is concerned. The net heating in the 
52^ 0 line is more extensive than that in the CO rotational 
lines. There is, of course, more energy available for absorp- 
tion in the planetary spectrum at 62 u than in the far infrared 
CO lines . 

For altitudes roughly 15 km above the altitudes of unit 
optical depth in 0 and CO the cooling given by (is) and (26) is 
less by more than an order of magnitude than that given by 
Bates' approximation. For several kilometers below the alti- 
tudes of unit optical depth there is an actual net cooling 
where the Bates' approximation gives no cooling. The form of 
the radiative emission terms, when expressed as the divergence 
of the radiative fluxes in the emission lines, therefore, 
deviates considerably from the Bates' approximation over a 
wide altitude interval. 
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A model atmosphere which is optically thin in 0 and CO 

has also been studied. The same temperature profile as that 

9 _3 

given by (3 2) was used, but a number density of 7 x 10 cm 
was assumed for both 0 and CO at 90 km. The net emission in 
this model was, at high altitudes, about 30% less than that 
given by the Bates' approximation, due to the absorption terms. 
At lower altitudes there was again a net heating. 

Summary and Conclusions 

Those terms in the thermal conduction equation which 
represent radiative emission by 0 and CO have been derived. 

The thermal conduction equation which includes these terms, 
as well as CO^ loss and ultraviolet heating, is the appropriate 
one to use in energy balance calculations for the thermospheres 
of the terrestrial planets. The difference between the equa- 
tion given in this paper and that which is usually employed in 
thermospheric investigations is the inclusion of infrared 

absorption terms in the former. The spontaneous radiative 

* 

emission at a given altitude is a function only of the local 
temperature at that altitude, but absorption depends upon the 
flux originating at all other levels and is thus a function of 
the run of temperatures throughout the thermosphere, as well 
as of the planetary effective temperature. The upper atmos- 
pheric energy balance is therefore expressed by an integro- 
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differential equation for the temperature profile T(z) in con- 
trast to the differential equation which expresses energy 
balance in the Bates' approximation. The net radiative emission 
as a function of altitude has been compared with that given by 
the Bates' approximation for a model Martian atmosphere. The 
magnitude of the deviation of the net emission terms from the 
Bates' approximation strongly suggests that the efficacy of 
this type of approximation be tested by including the absorption 
terms in thermal structure calculations for the thermospheres 
of the terrestrial planets. 
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Figure Caption 


Fig. 1. Comparison of the net emission, in ev/cm sec, 
in the 0 62 u line and the CO rotational lines with the radia- 
tive losses given by the Bates' approximation. Computed 
values of the various terms are based on a model Martian 
atmosphere which is optically thick in both 0 and CO. The 
curves are: (1) radiative heating and loss in O? (2) radia- 

tive heating and loss in CO; (3) radiative loss in 0 in Bates' 
approximation; (4) radiative loss in CO in Bates' approximation 
Curves (3) and (4) are terminated at the altitudes of unit 
optical depth* in O and CO. 
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